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g ; Abstract 

Self-assembled Plasmonic Nanoparticle Networks (PNN) composed of chains of 12-nm 
diameter crystalline gold nanoparticles exhibit a longitudinally coupled plasmon mode cen- 
tered at 700 nm. We have exploited this longitudinal absorption band to efficiently confine 
light fields and concentrate heat sources in the close vicinity of these plasmonic chain net- 
works. The mapping of the two phenomena on the same superstructures was performed by 
combining two-photon luminescence (TPL) and fluorescence polarization anisotropy (FPA) 
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imaging techniques. Besides the light and heat concentration, we show experimentally that 
the planar spatial distribution of optical field intensity can be simply modulated by controlling 
the linear polarization of the incident optical excitation. On the contrary, the heat production, 
which is obtained here by exciting the structures within the optically transparent window of 
biological tissues, is evenly spread over the entire PNN. This contrasts with the usual case 
of localized heating in continuous nanowires, thus opening opportunities for these networks 
in light-induced hyperthermia applications. Furthermore, we propose a unified theoretical 
framework to account for both the non-linear optical and thermal near-fields around PNN. 
The associated numerical simulations, based on a Green s function formalism, are in excellent 
agreement with the experimental images. This formalism therefore provides a versatile tool 
for the accurate engineering of optical and fhermodynamical properties of complex plasmonic 
colloidal architectures. 

Keywords : self-assembled plasmonics, two-photon luminescence, light-induced heating, gold 
nanoparticles. 

The optical properties of suspensions of individual colloidal metal nanoparticles have been 
scrutinized for more than a century, yet the rational design of their spectral properties by the tailor- 
ing of their size and shape has only become a topic of interest recently. 1 1 In the general context 
of plasmon-based technology, the next challenge for nanoparticle-based plasmonics is the control 
of the coupling between localized surface plasmons (LSP) by mastering the spatial organization 
of nanoparticles. 2 " 4 Self-assembly and templating principles have been proposed as the most ad- 
equate approaches to promote interparticle coupling and ordered multi- scale organization towards 
the design of specific nano-optical functionalities. 

In this context, we have recently reported the simple fabrication of complex and extended 
networks of interconnected chains of gold nanoparticles by a spontaneous one-pot self- assembly 
process driven by interparticle dipolar interactions. 2£& In suspension, these superstructures exhibit 
an overall globular size of typically 2-3 /im but individual chain segments are one-nanoparticle, i.e. 
12 nm, wide. As illustrated in[Q the optical absorption spectra of these superstructures, that we call 
Plasmonic Nanoparticle Networks (PNN), display not only a transverse plasmon mode (520 nm) 



but also a lower energy mode (700 nm) resulting from the strong coupling of nanoparticle surface 
plasmons between neighbouring nanoparticles along the chains and that we therefore called the 
PNN longitudinal mode 7 ' 8 . 

More recently, suspended PNN were deposited as intact and well-spread networks by deve- 
loping a specific substrate surface chemistry in order to control deposition and drying regime of 
the liquid colloidal solution. 9 

The combination of the micrometer-scale size of particle assembly and the 12-nm feature size 
of individual colloidal constituents confers to these self -assembled architectures the unique abil- 
ity to modulate very finely the spatial distribution of the electromagnetic field intensity. Once 
brought onto a solid surface, these colloidal superstructures could contribute to the improvement 
of applications as diverse as sensors, optical interconnects, enhanced nanoscale spectroscopy and 
microscopy. 

Concomitantly, the optical characterization down to sub-wavelength and molecular-scale fea- 
tures has made significant progress with the recent development of near-field 7,10-12 and single 
molecule microscopy, 13 i 14 electron energy loss spectroscopy (EELS) and imaging— or even suit- 
able far- field microscopy methods.-^"— In particular, two-photon luminescence (TPL), has been 
implemented in both near- field— ^ and far- field— ~— configurations. In the latter case, the TPL 
intensity can be related, in first approximation, to the square of the field intensity in the metal, i.e. 
the fourth power of the local electric field. 23 With such a fast-varying spatial dependency, far-field 
TPL microscopy reaches a spatial resolution as high as 200 nm. 24 Interestingly, optically-induced 
dissipation in plasmonic structures has also been characterized locally by far-field microscopy 
techniques such as Fluorescence Polarization Anisotropy imaging. 24 This indirect temperature 
measurement consists in measuring the loss of fluorescence anisotropy of a molecular probe in 
solution, upon linearly polarized excitation, when it crosses regions of elevated temperature due to 
heat dissipation in the vicinity of the illuminated plasmonic structure. When the entire plasmonic 
structure is illuminated while the FPA is probed locally, one measures the temperature distribution 
as detailed in references.— ^ On the contrary, when the illumination of the metal is confocal to the 



fluorescence excitation, then the FPA signal is related to the local heat source density. 
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Figure 1 : Typical extinction spectrum of plasmonic nanoparticle networks (PNN) revealing both 
transversal (520 nm) and longitudinal ( 700 nm) absorption bands. The arrow indicates the laser 
excitation of the longitudinal absorption band in our experiments. Inset: SEM image of a PNN. 
Scale bar: 500 nm. 



In this article, we report on the experimental observations and numerical simulations of near- 
field optical intensity and temperature mapping in the vicinity of gold PNNs. The large longi- 
tudinal band of these networks (FigfD is used to excite the coupled modes in the chains around 
730 nm with a hybrid microscope combining Two-Photon Luminescence (TPL) and Fluorescence 
Polarisation Anisotropy (FPA) measurements. 

The investigation of the plasmonic properties of lithographically designed metal structures by 
this combined technique was first introduced in reference 24 and the experimental set-up is illus- 
trated in|2]A Briefly, isometric gold nanoparticles were freshly prepared by the citrate reduction 



method 27 at a citrate: [AuCU] molar ratio of 5.2:1 and diluted to the required concentration with 
18 MQ. deionized water. The average diameter of the Au nanoparticles used in all experiments 
was 1 1.8 nm 1 1.1 nm. The assembly of the PNN was performed at room temperature by adding 
2-mercaptoethanol (MEA; HS(CH2)20H) to Au nanoparticle solutions at a Au nanoparticle:MEA 
molar ratio of l:5000.The nanoparticle chain assembly is characterized by a color change from 
pink to purple as the coupled mode emerges and was monitored by UV-vis spectrophotometry. 
The assembly is complete within 24 to 48 hours after mixing. Freshly prepared PNN are deposited 
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Figure 2: (Color online). (A) Schematic view of the combined TPL-FPA experimental set-up. (B) 
SEM image of a small PNN deposited on ITO/glass substrate. (C) TPL and (D) FPA images of the 
same PNN excited at 730 nm. 



on cleaned glass coverslips coated with 10 nm conductive Indium Tin Oxide (ITO) layer. 9 Well- 
spread networks were identified by FEG-SEM imaging (Figf2}B). To perform the TPL measure- 



ments, an IR Ti: Sapphire laser operated in pulsed mode is tightly focused on the sample (FigfJjB) 
with a high numerical aperture objective (NA=1 .3). The typical IR laser power is 1-3 mW and the 
pulse duration is 180 fs. In all the experiments, the IR laser wavelength was set to 730 nm. An 
image is obtained by raster scanning the sample in X and Y using a piezo stage and collecting the 
emitted light on an avalanche photodiode through a NA=1.3 objective (Figf2]C). 

To perform the temperature increase measurements, two laser beams are focused and over- 
lapped. 28 The IR TkSapphire laser is set in cw mode to locally heat the PNN, while a second lin- 
early polarized blue laser is used to excite the fluorescein molecules dispersed in a glycerol/water 
mixture surrounding the structure. The typical power of the blue and IR laser are 60 /iW and 4 mW 
respectively. The emitted molecular fluorescence is collected through the bottom objective and sent 
to two APDs via a polarizing cubic beam splitter to construct two fluorescence maps corresponding 
to 2 perpendicular emission polarization directions. The loss in fluorescence anisotropy is calcu- 
lated from the polarization difference signal and the associated temperature map is then obtained 
by using a FPA = f(T) calibration curve described elsewhere (FigfJt)). 28 To get the temperature 
increase as function of the IR laser location, the sample is raster scanned in X and Y as for the 
TPL measurements. For both types of experiments, particular care was taken to strictly avoid any 
damage of the nanoparticles by the laser-induced heating. The onset of melting could easily be 
identified in TPL by the rapid variation and eventually disappearance of the TPL signal within a 
few pixels. It usually occurred for nominal IR laser power of about 10 mW. The IR laser power 
used in all experiments was kept below 30 % of that needed for damaging the nanoparticles. 

Finally, the SEM images of the measured samples were resolved enough to extract positional 
maps of each individual 12-nm diameter nanoparticle from the studied PNN. These experimental 
maps were used to compute the simulated TPL and FPA maps. 

This theoretical analysis was based on a novel numerical tool involving the computation of 
temperature increase and TPL signal through a complete self -consistent scheme using both elec- 
trodynamical and thermal Green functions. 29 In the following, our model is first described in detail 
and experimental images are then compared to simulated ones for both TPL and FPA signals of 



one in several studied PNN. 

At the nanometer scale, the optical properties of small plasmonic objects deposited on a surface 
(small individual or self-assembled colloidal particles) can be investigated by adapting the coupled 
dipole approximation (CDA) method 30 to a planar geometry. 7 We first describe the unified formal- 
ism and numerical method that we developed in this context to compute the TPL and temperature 
distributions near PNNs, as depicted in|3] The network is composed of N identical subwavelength- 
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Figure 3: (color online) Schematic view of the experimental configuration in which a PNN made 
of 12-nm diameter nanoparticles is illuminated in normal incidence by a focus laser beam charac- 
terized by an electric field vector Eo(R w ,r,-, C0q), R w , r, and C0q being the beam waist position, the 
considered location for Eo and the incident field pulsation. 

sized nanoparticles (12 nm) located at the positions r,-. When the sample is illuminated by the 
electric field Eo(R w , r, coq) associated with a monochromatic laser beam of frequency 0)q centered 
and focused at R w , each nanometric component acquires an oscillating dipole moment driven by 
the local electric field E(R w ,r,-, (Do), which originates from the many-body interactions with the 
other particles: 



P(R w ,r f , coo) = a((0o) • E(R w ,r,-, (Oq) 



(1) 



where oc(coo) characterizes the dipolar polarizability of the metal particles. The many-body inter- 



actions between the nanoparticles can be introduced by writing N implicit linear equations: 

E(R w ,r ; -,coo) =Eo(R w ,r ; -,coo) 

+ J2S(r ; -,r ; -,(0 D )-a(oto)-E(R vv ,r ; ,(0 D ) , (2) 

i 

where S(r/,r/, C0q) represents the Green dyadic function of the bare substrate. 7 

The TPL signal arising from the gold nanoparticles is a non-coherent process characterized 
by a wide spectrum spanning from coo to 2(Do. This process is related to the intensity distribution 
|E(R vl .,r,-,ft)o)| 2 of local electric fields experienced by each particle. To solve this problem, we 
assume that the TPL signal can be modeled by summing the different intensities |E(R w ,r,-, coq) | 2 
rather than the electric field amplitudes E(R VV , r ; , (Oq) to account for the non-coherent nature of the 
TPL process: 
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N 
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(3) 



The overall square is intended to account for the two-photon absorption. 

In|H we present a sequence of theoretical (Figj4]\-B) and experimental (Figj4£-D) TPL images 
obtained by scanning the tightly focused Ti: sapphire laser spot with respect to a PNN composed 
of 1537 gold nanoparticles, the SEM image of which can be seen in the background of panels 
C and D. For each position of the laser beam waist R v , (incident wavelength Xq = 730 nm), the 
simulated TPL intensity has been computed from relation © by introducing the coordinates of 
the gold particle map plotted as overlay in black circles in panels A and B into the numerical 
code. The two incident polarizations, investigated along OX and OY directions, show an excellent 
agreement between experimental data and simulations, indicating that TPL measurements provide 
quantitative information on the light confinement near extremely small gold structures. Although 
the TPL signal seems to be more intense in regions of denser particle packing, this correlation 
does not account for the observed polarization dependency of the images. This polarization effect 
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Figure 4: (Color online) TPL mapping near a large PNN composed of 1537 gold nanoparticles 
(image size 2x2 /im). (A-B) Simulation images obtained from equation © overlaid with the 
nanoparticle map. (C-D) Experimental TPL images measured on this PNN overlaid with the SEM 
image of the PNN. Incident polarization is aligned along OX for (A) and (C) and along OY for (B) 
and (D). 



suggests that the local anisotropy of uniaxially coupled nanoparticles translates into anisotropic 
light confinement upon excitation of the PNN at 730 nm. 

Although limited by the intrinsic sensitivity and spatial resolution of our microscopy set-up, 
these observations illustrate how strongly coupled nanoparticles of small size but high crystallinity 
are able to laterally confine the electromagnetic intensity within a few tens of nanometers, while al- 
lowing its distribution along chains of a few hundreds of nanometers, as demonstrated theoretically 
elsewhere. 7 

Simultaneously, during the IR illumination, the PNN temperature rises because of Joule ef- 
fect induced inside the metal particles. The heat power Qi delivered by each nanoparticle can be 



computed from the electric field amplitude E;(R W , COq) using the relation: 24 

(o Im\a((Do)] . , .,, 

Qi(R w ,(Oo) = — \ K ^ i)i \E i (R w ,(O )\ 2 . (4) 



From this relation, the temperature rise AT(R W , C0q), expected at the laser location R w , can be 
deduced using the thermal Green's function: 29 

Ar M = ^f^>, (5) 

where K is the thermal conductivity of the surroundings. In order to take into account the influ- 
ence of the substrate conductivity, the image method has been used to compute the temperature 
distribution induced by the PNNs.— 

[5] shows a sequence of temperature maps resulting from the scanning of the same large PNN 
described in |4] through a 300 nm confocal spot of the near-IR and blue lasers . Although much 
less confined than the optical near-field intensity revealed in the TPL maps of HI the temperature 
distribution faithfully reproduces the PNN morphology. Noticeably, these thermal maps do not 
display any significant variation with the incident polarization (see maps [5fc and |5t)). In partic- 
ular, the three main branches of the network generate a same temperature rise because they are 
composed of many linear chains able to support the excited longitudinal mode irrespectively of 
the incident polarization. In addition, although we deal with extremely localized heat sources, the 
temperature distribution appears smoother than the near- field optical intensity map induced near 
the same system because of the thermal diffusion. 31 Two key differences between the generation of 
heat and optical TPL intensity should be emphasized. In the first case, as illustrated by Eq. ??, the 
temperature distribution generated at location R &$ by a point-like heat source located at R/, varies 
like the inverse of \R bs — Rft I while it varies like the inverse of the third power of this distance in 
the case of optical near-field. 

Furthermore, the TPL signal is sensitive to the square of the optical near-field intensity (cf. 
equation©, while the thermal process is just proportional to the global near-field intensity (cf Eq. 
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Figure 5: (Color online) Temperature mapping near the same large PNN shown in|4l (A-B) Simu- 
lation images obtained from ??. (C-D) Experimental images measured by the fluorescence polar- 
ization anisotropy of fluorescein molecules surrounding the PNN. Incident polarization is aligned 
along OX for (A) and (C) and along OY for (B) and (D). Calculated (A, B) and experimental (C, 
D) images are respectively overlaid with the positional network particle map and the SEM image 
of the PNN. 

??). This is clearly illustrated in both experimental and theoretical images of|H where we observe a 
significant variation of the TPL pattern inside the three branches of the PNN when the polarization 
is rotated by 90° in contrast to the invariance of thermal maps in|5] 

Interestingly, the excitation of longitudinally coupled modes in PNN involves electrically sep- 
arated nanoparticles in contrast to the continuous metallic stripes or nanowires prepared by litho- 
graphic approaches 24 or colloidal synthesis. As a consequence, the energy dissipation originates 
from each individual nanoparticle and is homogeneously distributed along the nanoparticles chains 
rather than localized where the excited nanorods plasmon mode induces a maximal current as 
shown in reference [24]. The self- assembly of the nanoparticles into globular micrometer-scale 
PNN therefore combines a number of specificities that could be advantageously exploited in the 
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context of hyperthermia cell treatment.— Although the full assessment of the hyperthermia effi- 
ciency of PNN on live bacteria, which is currently underway, lies beyond the scope of this article, 
the assets of our system already appear in the present results. Indeed, the heat production in PNN 
is obtained by addressing the longitudinally coupled mode centered at 700 nm that allows an ex- 
citation through the transparency window of biological tissues. Moreover, by choosing this wave- 
length, no individual stray particles or particle short chains (less than 4-5 nanoparticles) is excited, 
therefore obviating the parasitic heat dissipation away from the PNN. Moreover, the commensu- 
rability of the PNN with cell size, the biocompatibility of gold nanoparticles and the spontaneous 
affinity of PNN for the cell membrane (Figure SI in Supplementary Info) strongly suggest that 
optical energy could be dissipated directly at the PNN-wrapped cell walls without requiring an 
internalization step. 

In conclusion, this work demonstrates for the first time that large colloidal microstructures 
made of 12 nm diameter, crystalline nanoparticles self-assembled into single particle chains are 
able to effectively modulate the 2D spatial distribution of electromagnetic intensity, once deposited 
onto a substrate. Since the self-assembled networks are excited in the longitudinal coupled mode 
borne by the chains, these energy distributions can be simply tuned by changing the polarization 
of the incident light. The investigation of the near- field using TPL mapping provides an efficient 
means to characterize the intensity distribution at the scale of the entire superstructure but also 
reveals in-plane field modulation features which are only limited by the intrinsic spatial resolution 
of our TPL set-up, namely ca. 200 nm. Therefore, PNN are shown to have a unique potential for 
confining and enhancing the electromagnetic field down to sub- 100 nm length scale. 

Moreover, the measurement of the temperature distribution with our dual set-up has shown 
that heat is efficiently generated by exciting the longitudinal coupled mode and the generated heat 
is confined near the PNN. Understandably, the spatial distribution of the temperature elevation is 
more homogeneous and less sensitive to incident polarization than the TPL signal. Nevertheless, 
the optically-induced localized heat dissipation in our water-dispersible gold nanoparticle self- 
assembly holds promises in biomedical applications beyond the current capabilities of individual 
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colloids or substrate-bound metal structures. 

Finally, we have developed a model and simulation tool that accounts for most of the observed 
2D field map features. This is now a useful tool for the engineering of the optical near-field, 
for example for nanoscale optical information processing devices, using self-assembled structures 
made of crystalline metal colloids. The extension of our model and simulation tool to dissipation 
phenomena faithfully reproduced the observed temperature maps. 
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